Scaling Relations for Logarithmic Corrections 
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Multiplicative logarithmic corrections to scaling are frequently encountered in the critical behavior 
of certain statistical-mechanical systems. Here, a Lee- Yang zero approach is used to systematically 
analyse the exponents of such logarithms and to propose scaling relations between them. These 
proposed relations are then confronted with a variety of results from the literature. 
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Conventional leading scaling behavior at a second- 
order phase transition is described by power laws in the 
reduced temperature t and field h. With h — 0, the corre- 
lation length, specific heat, and susceptibility behave as 
£»(*) ~ I*! - ". C °o{t) ~ \t\- a , and X oo(t) ~ while 
the magnetization in the broken phase has moo(i) ~ |t| . 
Here the subscript indicates the extent of the system. At 
t = the magnetization scales (h) ~ h 1 ' 5 while 

the anomalous dimension r\ characterizes the correlation 
function at criticality. In the 1960's, it was shown that 
these six critical exponents are related via four scaling 
relations (see e.g. Ref. Q and references therein), which 
are now firmly established and fundamentally important 
in the theory of critical phenomena. With d representing 
the dimensionality of the system, the scaling relations are 
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In the conventional scaling scenario, © and © can, 
in fact, be deduced from the Widom scaling hypothesis 
that the Helmholtz free energy is a homogeneous func- 
tion 0. Widom scaling and the remaining two laws can, 
in turn, be derived from the Kadanoff block-spin con- 
struction Q and ultimately from Wilson's renormaliza- 
tion group (RG) Q . The relation 0) can also be derived 
from the hyperscaling hypothesis, namely, that the free 
energy behaves near criticality as the inverse correlation 
volume: foo(t) ~ , (i). Twice differentiating this rela- 
tion recovers Q. 

The scaling relations, 10) and J3J), were both rederived 
using an alternative route by Abe || and Suzuki |f| ex- 
ploiting the fact that the even and odd scaling fields can 
be linked by Lee- Yang zeros [jj . The locus of these zeros 
in the magnetic-field plane is controlled by the temper- 
ature. In the t > (disordered) phase this locus ter- 
minates at the Yang-Lee edge 0, the distance of which 
from the critical point is denoted by rvi,(t). At a con- 
ventional second-order phase transition ryi,(t) ~ t A for 



t > 0, and the gap exponent A is related to the other 
exponents through 



A = 
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Logarithmic corrections are characteristic of a number 
of marginal scenarios (see, e.g., Ref. and references 
therein). Hyperscaling fails at and above the upper crit- 
ical dimension d c and, while © holds there, it too fails 
above d Cl where mean- field behavior (which is indepen- 
dent of d) prevails. At d c itself, multiplicative logarith- 
mic corrections to scaling are manifest. Such corrections 
are found in marginal d < d c situations too 0,0. The 
g-state Potts model in d = 2 dimensions possesses a first- 
order transition for q > 4 and a second-order one when 
q < 4. The q — 4 case is also characterized by a transi- 
tion of second order, albeit with multiplicative logarith- 
mic corrections to scaling. Also in two dimensions, the 
Ising model with uncorrelated, quenched random-site or 
random-bond disorder offers another example of such cor- 
rections at a demarcation point. According to the Harris 
criterion [icl | , when the critical exponent a of the specific 
heat for a pure system is positive, random quenched dis- 
order is relevant (and exponents may change as disorder 
is added). If a is negative in the pure system, the critical 
behavior is not expected to be altered by such disorder. 
In the marginal case where a = 0, no Harris prediction 
can be made, and logarithmic corrections to the pure 
scaling behavior may ensue. These are some examples 
of the rich and disparate variety of systems displaying 
such phenomena and which have been hitherto studied 
individually. Given the ubiqiutous role that these loga- 
rithms play in such marginal cases, it is reasonable to ask 
if scaling relations for their exponents exist in analogy to 
0-0 above. 

Here three such relations, together with one for the 
Yang-Lee edge, are derived through the medium of par- 
tition function zeros and confronted with the literature. 
From the outset we mention that the scaling relations 
proposed herein do not all apply to the very special cir- 
cumstance of the Ising model in two dimensions and 
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its bond-disordered counterpart. These require special 
treatment beyond the general considerations presented 
here. For a scaling theory appropriate to the former, see 
Refs. an d references therein. With this in mind, we 

address the situation with the following scaling behavior: 
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while at t = 0, 



m oc (/i) - |ft| 1/5 |ln|ft| 



(11) 



In the thermodynamic limit the free energy may be 
written as 



/«,(*, ft) = 2Re f In [ft - ft(r, t)] ffoo (r, t)dr, (12) 



where i? is a cutoff, 3 M (r, t) is the density of zeros, with 
locus ft(r, i) = r exp [i0(r, t)]. If the Lee- Yang circle theo- 
rem holds the locus is given by = tt/2, # = 7T 0. While 
the validity of the Lee- Yang circle theorem is not assumed 
in what follows (it does not hold for the Potts model, for 
example), it is assumed that the small-i critical behavior 
is dominated by the zeros closest to the critical point, 
and that the locus of these zeros can be approximated 
by <j>(r,t) = (f>, a constant. 

The magnetic susceptibility is the second field deriva- 
tive of the free energy, and, at ft = Q [substituting 
r = xr Yh (t)] is 
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SooM) = Xoo(i)rYL(*)$ 
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up to additive corrections in ryi,{t)/R and where $ is 
an undetermined function of its argument. The ratio 
r YL (t) / R is sufficiently small near criticality so that these 
additive corrections may be dropped. Similar considera- 
tions yield for the magnetization 



moo ft) = Xoo(*)»"YL 
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Letting ft — > in l|15(l . and comparing to ©, recovers 
the scaling relation and yields 



A = /3- 7 . 



(17) 



Furthermore, fixing the argument of the function in 
|[TSJ gives t - ft 1 / A |ln|ft||- A / A from (0, so that itT5|) 
may be written 



moo^.ft) - ft 1 - 7/A |lnftp +7A/A * 



/ ft 
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Now taking t — > and comparing with (|ll|l recovers the 
known leading behavior for the edge JSJ, together with 
the correction relation A = S(S — 7)/((5 — 1). The former 
recovers ©, while the latter, with 1)170. gives 



/3(<5 - 1) = M - 7. 



(19) 



It is convenient at this point to define the cumulative 
distribution function of zeros as 



Goo (r,t) 



ryL(t) 



goo(s, t)ds = Xoo(t)r Yh (t)I 
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in which I(y) — <&{z)dz. Integrating \12\ by parts 
then gives the singular part of the free energy at ft = 0, 
/<»(*) = 2 J^ft) [Goo (r,t)/r] dr. Again substituting r = 
xrYL^), differentiating twice with respect to reduced 
temperature and comparing the resulting expression for 
the specific heat with JJJ) yields a = 2 + 7 — 2A and 
a — 7 + 2 A. From © and ©, the first of these is the 
scaling law (J2J. From l|17|l . the second can be conve- 
niently expressed as another relation between the correc- 
tion exponents, namely 



6 = 2/3-7. 



(21) 



Using these scaling relations, and fixing the ratio 
r/ryh(t) in l |2Qp and then taking the t — > limit, gives 
the critical cumulative distribution function to be 



Goo(r,0) ^ r (2-°0/A| lnr |fi-(2-a)A/z 
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Consider now a system of finite extent L, and let 
hj(L) = rj(L) exp (i<fij) be the j th zero there. The finite- 
size scaling (FSS) of first zero is expressible as 



ri(L) 
r Y L(t) 



&(Q) 



(23) 



in which (0) is the correlation length of the finite-size 
system at t = 0. On dimensional grounds, we may as- 
sume this quantity takes the generic form 



&(0)~L(lni)«, 



(24) 



having allowed for multiplicative logarithmic corrections. 
Recently, additional insights into the origin of FSS were 
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given in Ref. [12j. For a finite system, the cumulative 
density of zeros is simply the fractional number of zeros 
up to a given point, and we write 



Gl{t (L)) 
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(25) 



For large enough L, and at t = 0, this must coincide with 
the expression 122fl . In particular, it allows the scaling 
behavior of the lowest zero at thermodynamic criticality 
to be expressed as 



n(L) ~ir dA /( 2 -«) (InL) 



A-A5/(2-a) 



(26) 



Inserting iJBJ, l(TU|). J23J), and into J23) recovers Q 
and yields a new scaling relation for logarithmic correc- 
tions, namely 
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Hyperscaling corresponds to q = 0. Relations (|19fl and 
(I21|l but not l|27ll can be derived starting with a suitably 
modified phenomonological Widom ansatz [l3L fl4| . 
To summarize thus far, the three standard scaling laws 
(PJ have been recovered and three analogous relations 
for the logarithmic corrections (|19() . I|21(l . and (|27|l pre- 
sented. Furthermore, the standard formula 10 for the 
edge has been recovered and its logarithmic-correction 
counterpart is given in l|17|) . While the standard scal- 
ing laws for the leading critical exponents are well estab- 
lished, it is now necessary to confront the scaling relations 
for corrections with results from the literature, and a va- 
riety of models with logarithmic corrections arc examined 
on a case-by-case basis. 

The leading critical exponents for the 4-state Potts 
model in d = 2 dimensions were established in Ref. 0] 
as a = 2/3, /3 = 1/12, 7 = 7/6, S = 15, and v = 2/3, 
and their correction counterparts are [la . Il7| a = — 1, 
$ = -1/8, 7 = 3/4, 5 = -1/15, and v = 1/2. FSS of the 
thermodynamic functions are given in Refs. 0,0], from 
which q = 0. The standard scaling laws, of course, hold 
and one notes that the correction relations l|19|) . (|21|l . 
and (j23 hold too, while © and JHJ) give A = 5/4 and 
A = —7/8 for the edge. This latter prediction remains 
to be verified numerically. 

The upper critical dimension for O(N) symmetric <j)\ 
theories is d — d c — 4, where hyperscaling fails and the 
leading critical exponents take on their mean-field values, 
a = 0, /3 = 1/2, 7 = 1, S = 3, v = 1/2, A = 3/2. The 
RG predictions for the corrections are 0, 0, H3, 0] 
a = (4-N)/(N+8), $ = 3/(JV + 8), 7 = (N+2)/(N+8), 
6 = 1/3, v = (JV + 2)/2(/V + 8), A = (1 - N)/(N + 8), 
q = 1/4, and all of the correction relations (|17fl . (|19[) . 
(EH, and ((23 hold. 

The universality class of O(N) spin models can be ad- 
justed by introducing long-range interactions decaying 
(x being distance along the lattice), for which 



d c = 2a. The critical exponents for the TV-component 
system were calculated in Ref. [53 and are a = 0, 
/3 = 1/2, 7 = 1, 6 = 3, v = 1/<t, and obey the leading 
scaling relations. The Privman-Fisher form for the free 
energy was calculated in Ref. |23j, from which the RG 
predictions for the critical exponents could be verified 
and the logarithmic corrections observed. The logarith- 
mic exponents are a = (4 - N)/(N + 8), (3 = 3/(N + 8), 
7 = (N + 2)/{N + 8), <5 = 1/3, v = a(N + 2)/{N + 8). 
One observes that iflT)|) and lf2*T|) are obeyed, and lt2*7|) 
holds too if q — 1 /2a. This recovers the known value 
q = 1/4 for 0(N)(f)\ theory when a — 2, and leads to 
agreement with FSS in the long-range Ising case in two 
dimensions when a = 1 |24|. Furthermore, (JSJ) and i|17|) 
yield A = 3/2 and A = (1 - N)/{N + 8) for the edge, 
which are again identical to the 0(N)(f)\ values. The ex- 
pression 1)26(1 then gives that the first Lee- Yang zero of 
such a system should scale as ri(L) ~ L~ 3<T / 2 (lnL)~ 1//4 . 
This predicion for the Lee- Yang zeros of long-range sys- 
tems remains to be verified. 

Spin glasses, percolation, the Yang-Lee edge problem 
and lattice animals are all related to <fi 3 field theory. For 
each of these d c = 6 except the lattice animal problem 
which has d c = 8 ^25J. Ruiz-Lorenzo gave a compact 
description of the scaling of the correlation len gth, sus- 
ceptibility and specific heat for these models as [26j 
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The values of (a, 6) are (—4m, 1 — 3rn) for the m- 
component spin glass, (—1,-2) for percolation, and 
(—1,-1) for Yang-Lee singularities (which in d dimen- 
sions is closely related to the lattice animal problem in 
d + 2 dimensions). The mean-field values of the critical 
exponents for sp i n g lasses and percolation were calcu- 
lated in Refs. [271 l28j. respectively, as /3 = 1, 5 = 2, and, 
together with l|28|l . obey the usual scaling relations QJ- 
(0 . The correction exponents (|29|l satisfy the scaling re- 
lations lO and lEU provided that = 2(6 - a)/ (46 — a) 
and 8 = 6/(46 — a). In the percolation case these give 
(3 = 5 = 2/7, values which are in agreement with explicit 
calculations jll]. Also, © and (HZ} now yield A = 1 
and A = 2(6 - 2a)/(46 - a) for these models, while f2"7| 
gives q = 1/6 in each case. Ruiz-Lorenzo 's p rediction for 
this quantity is q = 1/3 26] while Ref. [23 contains an 
implicit assumption that q = 0. 

The strong universality hypothesis predicts that the 
quenched, disordered Ising model in d = 2 dimensions has 
the same leading critical exponents as in the pure case 
with logarithmic corrections to scaling EEl ] . In particular, 
Shalaev and later Shankar and Ludwig (SSL) gave [2II 
a = 0, j3 = 1/8, 7 = 7/4, S = 15, v = 1, a = 0, 
7 = 7/8, v = 1/2, with the specific heat predicted to be 
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double- logarithmically divergent |30|, and a more recent 
RG calculation gave $ = —1/16 and 5 = 0. 

Amongst the SSL values, that for 7 of the random- 
bond version has been the most clearly confirmed 
|33j ] . While the majority of published opinion favours 
the double- logarithmically divergent specific heat (see 
Refs. 0,13123 and references therein), there have been 
persistent claims in the literature that the specific heat, 
in fact, remains finite at criticality in the site-diluted 
model [3(| |32j • Compatability between the proposed 
scaling relations and the values (3 — —1/16, 7 = 7/8, 
5 = 0, and v = 1/2, is established if a = — 1 in this 
case. This value indeed leads to a finite specific heat 
in the random-site version and would neatly explain the 
persistent claims to that effect in the literature pjil |3tJ 
while still being consistent with the strong universality 
hypothesis. These values are also consistent, via (|27|) . 
with SSL's v = 1/2 provided q — 0, a value actu ally 
claimed on the basis of numerical evidence in Ref. |3J| 
(see also Ref. However, any value of v and q, dif- 

fering by 1/2, cannot be ruled out on the basis of (|27|l . 

Finally, the fact that the logarithmic divergence in 
the specific heat for the pure Ising model in d = 2 (see 
Refs. IHO]) does not directly fit into the scaling scheme 
proposed here is to do with special features of that model, 
which are shared by the random-bond version piU l32| . 
These special features are the vanishing of the specific- 
heat exponent a coupled with the property of self-duality 
and give rise to an extra logarithmic factor beyond those 
discussed herein. The apparent incompatability in these 
special cases with the scaling relations proposed herein is 
perhaps a reason why they have gone unnoticed as such 
before. 
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